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The general solution of the intertwining relations between a pair of 
Schrodinger Hamiltonians by the supercharges of third order in derivatives 
is obtained. The solution is expressed in terms of one arbitrary function. 
Some properties of the spectrum of the Hamiltonian are derived, and wave 
functions for three energy levels are constructed. This construction can be 
interpreted as addition of three new levels to the spectrum of partner poten- 
tial: a ground state and a pair of levels between successive excited states. 
Possible types of factorization of the third order supercharges are analysed, 
the connection with earlier known results is discussed. 

1. Introduction 



During last two decades the supersymmetrical (SUSY) method [1] for investigation of 
different problems in nonrelativistic Quantum Mechanics (QM) occupied its stable position 
among other more conventional approaches. In particular, it seems to be most adequate for 
study the properties of isospectrality (or almost exact isospectrality) of pairs of quantum 
models and to construct such pairs ("quantum design"). Correspondingly, the experience 
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of last years demonstrates that in the SUSY QM algebra the superinvariance of the Super- 
hamiltonian: 

[^,Ql = o (1) 

represents the main relation. Indeed, most of the known generalizations and applications 
of the standard Witten's SUSY QM [2] can be formulated as certain deformations of SUSY 
algebra, but Eq.(l), as a rule, remains unchanged (e.g. see [3], [4], [5]). In terms of the 
components of the Superhamiltonian H and supercharges : 

( H \ . , ( o\ 

; g+ = (g-)t= ; Q- = (g+)t, (2) 

\^ H j \Q- ) 

Eq.(l) takes the form of SUSY intertwining relations: 

HQ+ = Q+H; Q-H = HQ'. (3) 

Just these relations provide the isospectrahty (up to possible zero modes of Q^) of the 
partner Hamiltonians H, H and mutual connection between their wave functions: 

-^nix) = Q-^nix); ^n{x) = g+*„(x). (4) 

The first extensions of one-dimensionaP SUSY QM algebra, which preserve the inter- 
twining relations (3), concerned the supercharges of higher orders in derivatives and were 
called as the Higher order SUSY Quantum Mechanics (HSUSY QM), or alternatively, N-fold 
SUSY QM, or Non-linear SUSY QM [3], [4], [7], [8], [9]. 

Until now the intertwining relations were completely investigated for the cases of first and 
second orders. In particular, two kinds of second order supercharges were proven [4] to exist: 
reducible supercharges (arc equivalent to two successive first order supertransformations with 
real superpotentials) and irreducible supercharges (are not equivalent). 

Much less attention was paid in the literature to the transformations of third order 
in derivatives (see the papers [10], [11], [12]). In [10] the particular form of third order 

^On multi-dimensional extensions of Witten's SUSY QM see [6], [5] and references therein. 
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intertwining relations was investigated, which leads to a specific class of shape invariant 
potentials, while [12], [11] were mainly devoted to the construction of potentials with three 
lowest energy eigenvalues fixed. The main goal of this paper is to construct the general 
solution of intertwining relations (3) with supercharges of third order in derivatives. 

The paper is organized as follows. The most general solution of the intertwining relations 
(3) with supercharges of third order in derivatives is derived in Section 2. From this general 
solution of intertwining relations some particular properties of the spectrum are obtained 
in Section 3, the explicit expressions for three wave functions with fixed energy eigenvalues 
are written. Thereby one of the partner Hamiltonians, H, can be considered as quasi- 
exactly-solvable [13], [11]. The alternative "quantum design" interpretation is also given 
here: starting from H, involved in (3), one can build its partner H, whose spectrum includes 
three additional levels (a ground state ans a pair of levels between excited states). In Section 4 
the variety of possible factorizations of the third order intertwining operators is investigated. 
In Section 5 some particular cases of the general solution are studied, and its classical limit 
is discussed. 



2. The general solution of third order intertwining 

We consider the intertwining relations (3) on the real axis with the most general super- 
charges of third order in derivatives: 

Q+ = {Q-y^j:fn{x)d^; d = ^; xen (5) 

where without any loss of generality fsix) = 1 is chosen. Thus one obtains a system of nonlin- 
ear differential equations for real functions fni^), n = 0, 1, 2 and real potentials V{x), V{x) : 

V-V ^ 2/^, (6) 
/|-/2-2/i-3y = 3a, (7) 
2/2/1 - /r - 2/o - 3V" - 2/2^' = 0, (8) 



fi; + V"' + f2V" + f,V'-2fof^ = 0, (9) 

where a is an arbitrary real constant, and /' = df /dx. The derivative /q can be found from 
(8), and after insertion into (9) one obtains: 

(/i + Vy - 2h{f2 + Vy + 2f^{V - h)' + 4/o/^ = 0. (10) 

It is useful to introduce the function 

W{x) = V{x) + h{x) + a, (11) 

in terms of which other unknown functions in (6) - (9) will be expressed. Prom (7) and (11) 
functions fi{x) and V{x) can be written via W{x) and f2{x) : 

h{x) = m{x)-fi{x) + f^{x)- (12) 
V{x) = -2W{x) + fl{x)-f'^{x)-a. (13) 

In order to express everything in terms of one function W{x) it is useful to extract fQ{x) 
from (8) and (10). Then one obtains the equation: 

{w" + m"^)' - m'fj - 4/^ j f2W'dx = 0, 

and after its integration: 

W" + 6W'^ - 4/2 J f2W'dx = 2c, c = const. (14) 
Multiplying (14) by W and integrating again, one derives the relation between /2 and W: 
4( J f2W'dxf = + A{W^ -cW + d)] d^ const, (15) 
which gives the compact expression: 

W" + 6^-20 
= Tg ' 



4 



where g{x) is defined'^ by W{x) : 

g{x) = ±\yJW'^ + 4(1^3 _cW + d). (17) 

So, both the coefficient functions /„(x) of (5) for n = 0, 1,2 and potentials V{x)^ V{x) are 
expressed in terms of an unique arbitrary function W{x)^ thus giving the general solution of 
SUSY intertwining relations (3) of third order in derivatives by Eqs.(12), (13), (6), (8). 



3. Spectrum and wave functions 



Results of the previous Section lead to some consequences for the spectrum and wave 
functions of the Hamiltonian H. By direct (though rather cumbersome) calculations one can 
find that the operator Q~Q^ is represented as a third order polynomial of H with constant 
coefficients: 

Q-g+ = (// + a)^ - (// + a)c + d (18) 

Analogously, Q^Q~ is a third order polynomial of H with the same coefficients, indicating 
that the anticommutator of supercharges (2) is a third order polynomial of the Superhamil- 
tonian, i.e. one obtains the particular (third order) HSUSY QM [4] (see also [7] - [9]). 

It follows from the positivity of Q'^Q~ Eq.(18) that all eigenvalues of the Hamiltonian: 

H^,,{x) = En*„(x), (19) 

satisfy the inequality: 

{En + af - {En + a)c + d>Q. (20) 

We will now restrict ourselves by consideration of wave functions (19) which are simul- 
taneously the zero modes of the supercharge: 

Q+*fe(x) = 0. (21) 

"^We suppose that g{x) ^ 0. 
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We denote the energy eigenvalues, which correspond to these ^fe, by A^; A; = 0, 1, 2. It is 
clear that are the roots of the third order polynomial: 

{Xk + af - {Xk + a)c + d^O. (22) 

We notice that all three roots of (22) Aq < Ai < A2 are real if the constant parameters d and 
c satisfy the inequalities: 

c> 0; 27d^ < 4c^ (23) 

In this case (20) implies that Aq necessarily coincides with the ground state^: Aq = Eq, while 
Ai, A2 coincide with two neighbouring excited states: Ai = Ej; A2 = -E'j+i- Therefore all other 
(in general, an arbitrary number) energy eigenvalues En {n ^ 0, j, j + 1) may belong only to 
the following intervals: 

Ene{{Eo,Ej)U{Ej+,,oo)}. (24) 

Thus we deal with the so called quasi-cxactly-solvable [13] system, i.e. with known positions 
of the ground state energy and two additional (not necessarily first excited) energy eigenval- 
ues. In the case of only one real root Aq of (22), it realizes the ground state energy Aq = Eg, 
and other eigenvalues (if they exist) are to belong to the semiaxis E^ > Eq. 

These properties of the spectrum of H can also be interpreted from the "quantum de- 
sign" point of view. Let us suppose that the spectrum and eigenfunctions of the partner 
Hamiltonian H are known. Then if H is involved in the intertwining relations (3), one can 
construct the Hamiltonian H with the same (due to (4)) energy values and three (in the 
case of three roots above) additional levels Eq, Ej, Ej^i, one of which Eq becomes a ground 
state of H, and two other are inserted between a certain pair of excited states. In the 
case of one real root Aq only the ground state Aq = £^0 can be included in the spectrum of 
H. Of course, saying about new states ^„(x); n — 0,j,j + l, we have in mind that their 
wave functions are normalizable. All eigenfunctions of these added states will be constructed 
below, and their possible normalizability will be investigated. 

^We suppose that the corresponding wave function (see below) is normahzable. 
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By substitution of 9^ = {V{x) — H{x)) the Eq.(21) for A; = 0, 1, 2 is reduced to the hnear 
first order differential equation for ^n(^); — 0, j, j + 1. Its general solution reads: 

where En, n — 0,j,j + 1 is an arbitrary of the roots of (22). By means of (7) - (9) the 
nominator in the integrand can be rewritten as: 

V' + {V- En)f2 + fo^{W-En- a)/2 - W'/2 - g{x) 

to express the wave functions in terms of function W{x) : 

^n{x) =\ W{x) - {En + a) I V2 exp / (^-/, + ^J^^l + a) "^") ' ^^^^ 

where f2{x) and g{x) can be found from (16), (17). We remark that due to expression (26) 
one can keep the function W{x) fixed and study dependence of the wave functions \l/n(a^) on 
their spectral parameters En- Vice versa, one can fix the position of energy level, studying 
the class of potentials (varying the form of function W{x)) with the same values of En- 
It is appropriate to remind here that W{x) in (25), (26) is an arbitrary function, which 
must provide for the normalizability of eigenfunctions "ifnix)- To demonstrate that this re- 
striction for W{x) is not too strong, we give some examples of possible asymptotical be- 
haviour of W{x) at infinity and at finite singular point. 

1) Let W{x) has a growing power asymptotics at x ^ ±oo. Then a negative valued g{x) 
in (17) g{x) ~ —W'^^"^ and /2 ~ In this case "^{x) is normalizable if / W^^'^dx +oo 
for X ±oo. 

2) Let W{x) has an asymptotical x — > ±oo behaviour: W{x) — > 7 -|- ax~^, where 7(7^ — 
c) + d — 0- Then g{x) ~ ±^a{3^^ — c)x~^ and f2{x) ~ 2-\/37^ — ca'^x- Choosing a positive 
value for g{x), one obtains normalizable ^(x) ~ ex.p{—^/3^'^ — ca^^x^)- 

3) For d — and W{x) —a.x"'^, a > at x — > ±00 one has g{x) ~ i-y/cax"^ 
and f2{x) ~ ^^x l2y/oL- For the negative sign the wave functions are normalizable: ^ ~ 

exp ( — y^X^ /4 y^) . 
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4) If W{x) has a singularity at x — 0: W{x) — > ax"'^, then g{x) ~ ±a^/a + lx~^, and 
/2 ~ ±3y/a + Ix"^ /2. For the negative sign ^ x j ^. ^(x) is normahzable at 
X = for a > 0. 

5) For a pole type singularity of W{x) ~ ax^^ at a; ^ one can take the negative sign 
of g{x) ~ —ax~'^/2 and f2{x) ~ —x^^ to obtain normahzable ^(x). 

In all cases it is necessary to watch nodes of the function g{x), which generate the 
singularities of f2{x), in order to keep the normalizability of wave functions under the control. 

For illustration of the approach we consider a particular example with W{x) =7 + j;^, 
where, according to item 2) above, we must choose a{3^'^ — c) > and d — 7(0 — 7^). Let us 
take in addition 7 = —a{a + 4)/12 to ensure that the function g{x) is an entire function: 

. . a'^x(x'^ + — ) ^ . . ax a - 2 2x 



2(l + a;2)2 ' ■'^^ ^ 4 Ax 1 + x^' 
Then the partner potentials have (up to a common additive constant) the form: 



~ , , o (a^ - 4) 1 4 8 a , , 

The singular harmonic oscillator potential (27) on the entire axis-^ was studied in the liter- 
ature (e.g. see [14]). We would not like to make here the detailed study of this potential 
our aim, pointing out the paper [15], where this kind of singular potentials was investigated 
just in the frameworks of SUSY QM and unbroken shape invariance. Summarizing [15], the 
spectrum for the transition (see for definitions [16]) "soft" (i.e. for — | < ^""^^g""^'* < |) 
potential (27) was shown to consist of two equidistant sequences of harmonic-oscillator-like 
levels, and corresponding wave functions were also constructed. 

The condition of "softness" of both transition potentials (27) and (28) means that a e 
(0,4). Thus, the third order intertwining of V{x) and V{x) ensures that the spectrum and 
■^Though over this paper we deal with the problems on a whole line, it can be directly restricted to a half 
axis a; > by imposing suitable conditions for behaviour of wave functions at the origin. 
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wave functions^ of the potential (28) coincide (as usual, up to zero modes of Q^) with those 
of the partner potential (27). 

4. The factorization of intertwining operators 

The important question concerns the possible factorization of the general intertwining 
operator (5) onto the lower order multipliers. We start from the factorization onto three 
operators: 

Q+ = -qtqtqt, (29) 
where q^; i = 1,2,3 are of the first order in derivatives'^: 

q+ = -d + Wi{x). (30) 

Eq.(29) leads to the following system of nonhnear differential equations between superpo- 
tentials Wi{x) and coefficient functions fn{x) : 

/2 = -iWi + W2 + Ws); (31) 

/i = -{W2 + W3y + Wi{W2 + W3) + W2W3-w;,; (32) 

/o = {W2W3y -WiW2W3-W^ + WiW^. (33) 

One can take W2 + W3 from Eq.(31), insert it into Eq.(32) and use the result in Eq.(33). In 
such a way one obtains the nonlinear second order differential equation for Wi{x) : 

{W^ - W[y - {W^ - W[){W, + f2) + Wi{2f^ - /i) - + /( - /o = 0. (34) 

Eq.(34) can be linearized by defining Wi{x) as a logarithmic derivative of a function y{x): 

y"'-f2y" + {fi-2f!,)y' + {f[-f^-fo)y^0; W^i(x) ^ (35) 

^The nonsingular (a = 2) particular case of (28) was discussed in [11]. 

^We choose the signs in operators the same as in the conventional first order SUSY QM [1] . 
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Thus an arbitrary of solutions of (35) gives the real superpotential Wi in terms of real 
functions fn{x). Then one can use (32), (33) to find also W2 and W2, in terms of /„. 

It is necessary to remark that for the real valued solution Wi one can consider two 
different opportunities for W2 and : 

1) real valued W21 W^] and 

2) mutually conjugated complex W2 = with constant imaginary part*. 

The first option can be refered to as realizing a maximally reducible supercharge (su- 
pertransformation) . The first order multipliers create a chain ("dressing chain" [17]) of 
Hamiltonians (—9^ + {W^ + e^), such that the neighbouring Hamiltonians (up to a con- 
stant) in the chain are intertwined by q^. This procedure is a particular case of the familiar 
construction of reducible Higher order SUSY QM [3], [4], [11] by gluing. In this case the 
polynomial (22) has three real roots En, n = 0, 1, 2. 

The second option corresponds to partially reducible supercharges, where = qiM2s 
with irreducible second order operator with real coefficient functions: 

M23 = qtqt = {-d + W2{x)){-d + W*{x)). (36) 

Though the superpotentials W2, are chosen here to be complex, the coefficient functions 

/„ and potentials V ,V still are kept reaP. Just this kind of intertwining operators was 
investigated in [3], [4]. In this case (22) has one real and two mutually conjugated complex 
roots Ajfc. 

'In this case the third order polynomial (22) has one real and two mutually conjugated roots. 
■'The complex valued /„ (with real V,V) would lead to additional intertwining relations of lower order 
(see for details the preprint in [10] and [9]. 
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5. Some particular cases of intertwining 



5.1. Quasi-exactly-solvable models with three lowest levels given 

The particular case of the first (maximally reducible) option above has to be compared 
with the construction of the paper [11], where the intertwining operators of third order 
were useed to build the quasi-exactly-solvable Hamiltonians, with three lowest values of 
bound state energy given. In [11] from the very beginning the intertwining operator had 
the factorized form with real superpotentials Wi{x). Those superpotentials were generated 
by the groung state wave functions '^o\x) of three consecutive Hamiltonians H^'^^; i — 
1, 2, 3, included in the chain [4] of three supersymmetric transformations of first order with 
intermediate gluing (up to a constant shift): 

i/«*«(x) = £;^^*?(x); ^S\x) = exp(- 1 W,{x)dx). (37) 

It was shown in [11] that all three superpotentials Wi{x) and corresponding potentials can 
be expressed in terms of one function U {x) with definite strong limitations in its analytical 
properties originated from the restrictions on logarithmic derivative of normalizable ground 
state wave functions (37): JWi{x)dx has no singularities, and asymptotically sign Wj(x) = 
±1 at X — > ±oo. It can be shown by rather long calculations that the function U{x) of [11] 
and W{x) in Section 2 of the present paper coincide up to the constant if the condition (23) 
is fulfilled. Thus the solutions described in [11] can be obtained as particular cases of the 
general construction presented above, if one will restrict himself by suitable strong conditions 
for analytical and topological properties of functions and for relations between constants. In 
the general construction of the present paper both the singularities and the complexity of 
functions Wj are allowed if they preserve potentials V, V to be real and not too singular. 

We remark that each of two kinds of factorization of depends on one arbitrary (up to 
some discrete restrictions) function. It is obvious that the alternative to option 2) partially 
reducible factorization of Q'^ exists, if one chooses W3 real and Wi = W2 with constant 
imaginary part. The third similar opportunity with W2 real and Wi = is equivalent to 
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the previous ones due to commutation [qi, — [q^, qt] — 0. 



5.2. Models with shape invariance of third order 

We will now consider the particular case of intertwining relations (3) of third order where 
H and H coincide up to a constant 2A : 

HQ+ ^Q+{H + 2X). (38) 

This (shape invariant) sort of intertwining appeared in [10] in a natural way when a pair of 
Hamiltonians H, H participated simultaneously in two intertwining relations - of first (by 

operators q^) and second (by operators M~) order in derivatives. 

q+ ^ -d + u;(x)); ^ d'^ + 2(j){x)d + b{x) + 2(l)'{x)), (39) 

where (see [10], [4] for details) 

= ^^^ + 6^'(rc) + 8WW+2(AV-(A + a))^(x) + ^ (41) 

The specific shape invariance [18] property of the Hamiltonian H, realized by (38), allowed 
to find its spectrum algebraically [10] by analysis of zero modes of and Q~ : 

Q+ = q+M-; Q" = M+q'. (42) 

Comparing coefficients of d'' in (42), one obtains the relations: 

f2{x) = 2(j){x) -ij{x); 
fi(x) = 40'(x) - 2a;(x)0(x) + b(x). (43) 

The function f2{x) is fixed by (6) and (38): 

f2{x) = —Xx; u!{x) = 2(f>{x) + Xx. 
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It follows from Eqs.(12), (14) and (43) that 



W{x) = (p'{x) - 24>^{x) - 2Xx(i){x) + 



2X + a 
3 



and, due to (14), has to satisfy nonlinear differential equation: 




(44) 



Indeed, it can be shown by direct calculation (using (41)) that (44) is satisfied for c = 



p. Thus the shape invariant case (38) is included in the general solution of Section 2. 



5.3. Models with intertwining of first and third orders 

The third order intertwining of the previous Subsection was based on the double in- 
tertwining of the same pair of Hamiltonians H,H - hj operators 5+ of the first and M" 
of the second orders. Now we will consider the double intertwining by operators q'^ of first 
order and Q''" of third order in derivatives. The first intertwining gives: 

H = q+q- = -d'^ + u'^{x) -uj'{x) + e; H = -d'^ + uj^{x) + uj'{x) + e; e = const, 

while the second one (see (6), (13)): 

V{x) = -2W{x) + fl{x) - f^{x) -a = V{x) - 2f^{x). 



Two cases have to be considered separately: 7 = and 7 7^ 0. 

For 7 = from (45) W{x) — —\{a -\- e) = hi — const. Then from Eqs.(6) - (10) one 
can derive that 



Thus, 



uj{x) = -f2{x) --f; = const; 
27/2(2;) = -2W{x) - (a + e + 7^). 



(45) 



fo{x) = SKf2{x) 



h{x)^?,K-fl{x)+f^{x)- 

V{x) = f^ix)-f!,ix)-2n-a. 
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Then the operator Q'^ is reduced ("modulo Hamiltonian if") to the first order operator 




(k - a - H)q^. 



More interesting case is for 7 0, when 

f2{x) = — 7~^iy(a;) — v = const. 

Substitution of this expression into Eq.(15) gives equation, which includes the function W{x) 
only: 



In general case its solution can be expressed in terms of elliptic functions, but some particular 
values of constants (e.g. v — c — d — ^)v^ (46) give simpler solutions as well. 

It can be shown that now the operator is expressed (again "modulo Hamiltonian if") 
in terms of first order and second order operators 



Therefore the same pair of Hamiltonians is intertwined both by the first order and the 
second order operators. This result is in accordance with the general statements of the 
recent paper [9] concerning "the optimal set of two basic SUSY generators of even and odd 
order". One can check also from the intertwining relations that the third order operator 

R = (L^ Yq^ commutes with H, playing the role of its symmetry operator (see discussion in 
the preprint version of [10] and [9]). 

5.4. The classical limit 

It is known that sometimes results obtained for the quantum (supersymmetrical) systems 
can be successfully used to investigate its classical prototype (e.g. see [19] on two-dimensional 
second order SUSY Classical Mechanics). We will consider here the consequences of third 
order SUSY QM for the classical limit h ^ 0. With this object we have to restore the h in 



(W'f = -f-^W^ + 4(1/7-^ - 1)W^ + Aiy^W^ + AcW - Ad. 



(46) 



g+ = L+ - Hq+; L+ 



7-^92 + {W- a)d -Vu^ J'. 
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all formulas of Section 2. It can be done by rewriting differential operators in terms of the 
momentum p ~ —ihd. The mnemonic rule is to multiply each derivative by the multiplier h 
in Eqs.(5) - (10). Therefore, in the classical limit ^ — > one obtains from (6) - (9): 



V-V = 



{fl-2h-W)' = 0; 
fUi-fo-f2V' = 0; 
fiV'-2fof!, = 0. 



(47) 



The general solution of this "classical system" of equations can be obtained from its quantum 
counterpart of Section 2 by mnemonic rules formulated above. It is useful to compare this 
problem with discussion in the paper [8], where the function B{x, ip) — J2 bs-niw)"' plays the 
role of classical analogue of the quantum operator (5+ with &3_„(x) = fn{x). These coefficient 
functions have to satisfy the differential equation (Eq.(2.24) of [8]): 



which was not solved in the framework of pure classical theory (see details in [8]). However 
Eq.(48) coincides exactly (up to some normalization factors) with Eq.(47), whose general 
solution can be found analytically from the quantum solution of Section 2 by ^ — > hmit. 
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b'^{x) + (5 - n)br,-2{x)V'{x) - L{x)hn-i{x) = 0, 



(48) 
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